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Abstract. In the LHC collider a significant rate of events with double parton 
scattering is expected. This will be the leading mechanism for production of two 
bb pairs. We estimate the probability of binding two b quarks into a diquark and 
the probability of dressing this diquark into a bbud (IS P = 01 + ) tetraquark. 
Calculations shows that that this bound state of two B mesons is stable against 
the strong interaction and has a life time of the order of ps. We estimate that 
the production rate at luminosity L = 0.1 events per second will be about 6 
tetraquarks per hour or more. 

1 Introduction 

In the new large hadron collider at CERN on can expect a sizable rate of events 
where two heavy b quarks and two heavy b antiquarks are simultaneously pro- 
duced. This makes it possible that two heavy quarks or antiquarks combine to 
create a double heavy hadron after they dress with light quarks. 
There are two mechanisms for creation of two heavy quark- antiquark pairs. 
The first is by a single parton scattering where initially only one heavy quark- 
antiquark pair is produced and the second pair is then created by the frag- 
mentation of one of initial quarks. The second mechanism is a double parton 
scattering. Calculations show that the double scattering although power 
suppressed gives at the center of the mass energy of 14 TeV a four times larger 
cross section for the production of two b-quark pairs. For the creation of a 
heavy diquark one requires that the two quarks are close together in momen- 
tum space. Since in the fragmentation the two b quarks tend to fly in opposite 
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directions this also confirms that the dominant mechanism would be two parton 
scattering. The creation of the double heavy hadron is then a two step process. 
First two heavy quarks which are close together in momentum space bind into 
a diquark. We can estimate that the probability for this is proportional to the 
overlap of the wave functions of the quarks which are considered to be free 
with the wave function of the diquark The second step is then dressing of this 
diquark with the light quark into a doubly heavy baryon or dressing with the 
two light quarks into a heavy dimeson. 



2 Phenomenological estimate 

Let me first give a very simple phenomenological estimate why a BB dimeson 
should exist with a binding energy of about 100 MeV and why a similar DD 
dimeson is unstable against strong interaction 2 . In this estimation we consider 
the heavy dimeson to be a bound state of two heavy b-quarks combined in a 
diquark with the spin 1 and two light antiquarks with the quantum numbers 
spin 0, isospin 0. The starting point is that the light antiquarks around the 
heavy diquark behave as in A^. This would be exactly true in the limit where 
the mass of the b quark goes to infinity and the heavy diquark is point-like. This 
correspond to assumption that we can neglect the size of the heavy diquark in 
dimeson. The second assumption is that the interaction between heavy b quarks 
in a colour antitriplet state is half as strong as the interaction between heavy 
quark and antiquark in heavy meson. With this assumption we can estimate 
the binding energy of heavy dimeson. Since the spectroscopy of heavy mesons 
is quite well described with nonrelativistic potential models we see that the 
binding energy for a system with half as strong interaction Vqq is half the 
binding energy for a system with original interaction Vqq and twice lighter 
particles (toq — > jrq/2). 
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We obtain E(itiq/2) by plotting the binding energy of heavy mesons as a 
function of the reduced mass of the system j2J . We do that for different choices 
of constituent quark masses, to estimate the uncertainty due to quark masses. 
The plot is very smooth and we estimate the binding energy of the heavy bb 
diquark E(rrib/2)/2 to be —390 ± 15 MeV. If we combine this with the masses 
of heavy mesons and Ab baryon we can obtain the phenomenological estimate 
for binding energy of the heavy dimeson ZlTbb 

AT hh = m Ab + [m T - E{m h ) + E(m b /2)] /2 - m B - m B « = -130 ± 15 MeV 



while the heavy dimeson with the two c quarks is unbound. This agrees well 
with some previous four body calculations [3] 0] in the constituent quark model. 
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3 b-quark pair production 

Since at LHC there will be a large parton luminosity we can hope that the 
production rate of double heavy dimesons will be high enough to make detection 
possible. The dominant contribution to b-quark pairs production will be double 
parton scattering - the gluon fusion (gg — ► bb) 2 . To calculate this cross section 
one needs to make an assumption about the two body distribution function r 

r(x 1 ,x 2 ,d) = G(x 1 )G(x 2 )F(d), (1) 

where G(x) is the one-body parton function while the parton pair density is 
described by the normalized function F(d). The distance of the two partons be- 
longing to the same proton is d and the x\ and xi are their fractional momenta. 
Doing so, we neglect the correlations in momenta. Then the cross section <jd 
for the double gluon fusion has a form 

^(bbbb) = SiW! a -$ = f d HF{bf, (2) 

where a s (bb) is cross section for a fusion of two gluons into a b quark-antiquark 
pair. For a e ff we take the value measured in the experimental study of double 
parton collision by CDF at Fermilab in the reaction pp — » j/ir° + 3 jets + X 
This experimental value is 14.5 mb. Doing this we assume that also some 
correlation between quarks in the same proton are taken into account which 
were neglected in Q . There is also same uncertainty about the single scattering 
cross section a s since the higher order corrections are large and they enter 
squared into the cross section for double scattering (J2J). For details see £Q. 
The cross section calculated then from (|2J), taking into account that only the 
pseudorapidity region 1.8 < rj < 4.9 will be covered by forward LHCb detector, 
is 85 nb. For two b quarks to combine into a diquark it is very important 
that they are close together in momentum space. Calculations gives that the 
cross section is approximately proportional to momentum volume da/d 3 p « 
0.4nb/GeV 3 up to p < 2GeV. 



4 Formation of the diquark and dressing 

We assume that the production of the two heavy b quarks in double parton 
scattering is uniform inside the protons. Furthermore we take an instantaneous 
approximation so that the amplitude M(p) for the formation of the diquark is 
simply given by the overlap M(p) = (4>(p)\ip) of the wave function of the two 
b quarks <j>(p) with the wave function of the diquark ip. This wave function is 
calculated in constituent quark model used also in [3] and 0] . So the production 
cross section is given by 

*= [ d 3 p^-M 2 {p) = 0.15 nb (3) 
J d^p 



Dressing of a heavy diquark is governed by the quantum chromodynamics. But 
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since fragmentation of diquark into a hadron is impossible to calculate using 
perturbative QCD methods, we estimate the probability from experimental 
data obtained at Fermilab and at LEP experiment 6 . Their results show that 
the probability that the b quark will be dressed with u or d antiquark is 0.38 ± 

0. 02. with s antiquark 0.16 ± 0.03 and the probability for dressing with light 
quarks to form At is 0.10 ± 0.03. In our estimation we compared the heavy 
dimeson with the A\> baryon where the role of diquark in the first case is 
played by the heavy b quark in the second. So we assume that the diquark will 
get dressed into a dimeson with the same probability as heavy quark is dressed 
into a lambda b baryon - that is 0.10 ± 0.03. The cross section for the heavy 
dimeson production is thus 

a = 0.015 nb 

which will give at the expected luminosity the production rate of 5-6 events 
per hour. 

5 Decay and detection 

We expect that the bbud tetraquark will be stable against strong and electro- 
magnetic decay and will decay only weakly. The main channel would be the 
decay of one of the b quarks into c quark followed by the weak decay of the 
second b quark. But this channel will be difficult to distinguish from the chan- 
nel of two unbound B mesons which will be the main background contribution. 
The more characteristic decay channel would be the the formation of the T and 
light meson through b — > b oscillation inspired by B° — ► B° oscillations. But 
this oscillations are here negligible since the tetraquark is not degenerate with 
final state of two unbound mesons. We are still searching for new ideas for the 
detection of doubly b tetraquarks. 

References 

1. A. Del Fabbro and D. Treleani, Phys. Rev. D 63, 057901 (2001) 

2. D. Jane, M. Rosina: Few-Body Systems 31, 1 (2001) 

3. B. Silvestre-Brac, C. Semay: Z. Phys. C57, 273 (1993) 

4. D. M. Brink, Fl. Stancu: Phys. Rev D57, 6778 (1998) 

5. F. Abe et al. (CDF Collaboration), Phys. Rev. D 56, 3811 (1997) 



6. T. Affolder et al. (CDF Collaboration), Phys. Rev. Lett. 84, 1663 (2000) 



